A near-infrared study of the star forming region RCW 34 
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ABSTRACT 

We report the results of a near-infrared imaging study of a 7.8 x 7.8 arcmin 2 region centered 
on the 6.7 GHz methanol maser associated with the RCW 34 star forming region using the 
1.4m IRSF telescope at Sutherland. A total of 1283 objects were detected simultaneously in J, 
H, and K for an exposure time of 10800 seconds. The J-H, H-K two-colour diagram revealed a 
strong concentration of more than 700 objects with colours similar to what is expected of reddened 
classical T Tauri stars. The distribution of the objects on the K vs J-K colour-magnitude diagram 
is also suggestive that a significant fraction of the 1283 objects is lower mass pre-main sequence 
stars. We also present the luminosity function for the subset of about 700 pre-main sequence stars 
and show that it suggests ongoing star formation activity for about 10 7 years. An examination 
of the spatial distribution of the pre-main sequence stars shows that the fainter (older) part of 
the population is more dispersed over the observed region and the brighter (younger) subset is 
more concentrated around the position of the 08.5V star. This suggests that the physical effects 
of the 08.5V star and the two early B-type stars on the remainder of the cloud out of which 
they formed, could have played a role in the onset of the more recent episode of star formation 
in RCW 34. 

Subject headings: stars: formation — stars: pre-main sequence — ISM: individual (RCW 34) 



1. Introduction 

Studies on star formation can very broadly be 
classified as either the study of the process of star 
formation as such or as the study of the larger scale 
processes that trigger star formation in molecu- 
lar clouds and which determines the mode of star 
formation, ie. clustered or distributed star for- 
mation. Although large strides have been made 
over the last couple of decades in understanding 
the process of star formation, an understanding 
of what determines the mode of star formation, 
which is equally important to understanding the 
process of the formation of individual stars, is still 
lacking. The two aspects are, obviously, not in- 
dependent of each other and together play a role 
in determining important properties such as, for 
example, the stellar initial mass function and the 
star formation efficiency. To achieve the goal of a 
theoretical understanding of the physical process 
determining the mode of star formation, empirical 
information on a sufficiently large sample of em- 



bedded young clusters is necessary ( Ladal 2010h . 
The present near-infrared study of RCW 34 at- 
temps to add to this effort. 

RCW 34 is a southern high mass star form- 
ing region that has been studied at irregular in- 
tervals over a number of decades. The number 
of observational studies that was directly aimed 
at investigating RCW 34 is quite small. We re- 
view those that are relevant to the present study 



RCW 34 was first catalogued by iRodgers et al 
( I960! ) as one of the southern Ha emission re- 
gions with a diameter less than 2x2 arcmin- 
utes. The exciting star of RCW 34 also appears in 
the catal ogue of stars in reflectio n nebu lae com- 
piled by Ivan den Bergh fc Herbst (1975). Based 
on this catalogue the exciting star in RCW 34 
is also refered to in the literature as either Vela 
R2 25a or vdBH 25a. Following up on t he cata- 
logue o f Ivan den Bergh &: Herbsd (|1975l ). iHerbst 
(|l975af ) did UBV and MK spectroscopy on the 
stars in the catalogue and classified the exciting 
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star in RCW 34 as of spectral type 09. In a fur- 
ther analysis of the UBV data, iHerbstl (|l975bl ) did 
main-sequence fitting on the stars of Vela R2 as- 
sociation. This resulted in a distance of about 870 
pc to the association. However, the exciting star 
of RCW 3 4 was found to l ie well below the ZAMS, 
which l ead IHerbstl (Il975b ) to consider it to be pe- 
culiar. iHerbstJ (jl975t¥ explained the position of 
RCW 34 well below the ZAMS as due to it being 
either an O star located further than the Vela R2 
association or an O-type pre-main sequence star 
with a " grey" circum stellar shell associated with 



also ide ntified ind i catin 
tivity. iBik et all jioT 



Vela R2. Vittone et al 



9871) w a s the fi rst to fol- 



Herbsd (|1975b[ ) that the 



low up on the finding of 
exciting star of RCW 34 lie well below the ZAMS 
for the Vela R2 association. These authors also 
concluded that vdBH 25a is an 09 type star but 
that the distance to the star is uncertain. 

The first really comprehensive study of RCW 
34 aimed to answer the questions on the na- 
ture of vdBH 25a and its distance is that of 
Hevdari-Malaveri (1988). Using a v ariety of obser- 
vation s and data from the literature, Hevdari-Malaveril 
(|l988l ) concluded that RCW 34 is located at 2.9 
kpc, well beyond the Vela R2 association, and 
that its exciting star is a heavily reddened (A„ 
= 4.2) 08.5V ZAMS star. There is therefore 
nothi ng peculiar about the exciting star of RCW 
34. Hevdari-Malaveri (|l988l ) also interestingly re- 
marks that the presence of an H2O maser might 
signify ongoing star formation in RCW 34. 

The most recen t study of RCW 34 is that 
~ f)2010h . 



of Bik et al 



These authors made near- 
infrared H- and K-band observations on a 104" x 
60" field centered on the 08.5V ionizing star, us- 
ing the Integral Field instrument SINFONI on 
UT4 of the VLT at Paranal, Chile. In addition 
they also used arc hival data obtained with IRAC 
(|Fazio et al.1 12004) on boa r d the Spitzer satellite. 
Using this data IBik et "all (|2010l ) identified three 
regions associated with RCW 34: The first is a 
large bubble detected in the IRAC images in which 
a cluster of intermediate- and low-mass class II 
objects is found. The second is the H II region 
located at the northern edge of the bubble and 
which is ionized by three OB stars. The third re- 
gion is a photon-dominated region north of the 
H II region and which marks the edge of a dense 
molecular cloud traced by H2 emission. Several 
class 0/1 objects associated with this cloud were 



recent star formation ac- 
revised the distance to 
RCW 34 to 2.5 ± 0.2 kpc and derived an age es- 
timate of 2 ± 1 Myr from the properties of the 
pre-main sequence stars inside the H II region. 

The observations presented here differ from 
earlier near-infrared (N IR) observations , and in 
particular from that of IBik et aL ( 2010h in that 
we imaged a 7.8' x 7.8' region centered on the 
6.7 GHz methanol maser located at coordinates 
a(2000) = 08 h 56 m 24 s .9 and 

5(2000) = -43°05'40".4. This is significantly 
larger than the 104" x 60" region covered by 
IBik et all (|2010l ). This much larger field gives us 
the opportunity to also investigate possible star 
formation activity especially to the north of the 
shocked region where the H II region interacts 
with the molecular cloud and crea ted a photo- 
dissociation region (Bik et al.l [ioToh . Except for 
the class II methanol maser which indicates the 
presence of a very young high mass star there ap- 
pears to be no other signs of current high mass 
star formation activity to the north of the H II re- 
gion. This does not exclude the presence of lower 
mass stars, however. 

Within the framework of a broader view on star 
formation as outlined above, our results should be 
seen a s complimentary to that of eg. 



Bik et al 



(2010) and others which can help to understand 
the star formation history of RCW 34. Our aim 
here is to present some basic properties of the em- 
bedded cluster associated with RCW 34 as derived 
from our NIR observations. 

2. Observations, data reduction, and cali- 
bration 

The observations used in this study can be re- 
garded as archival data in the sense that the data 
were not collected by the authors themselves for 
the primary purpose of studying RCW 34 as is 
done here. The data used here was collected by 
GG Nyambuya in 2005, however, for a completely 
different purpose. This caused the data to be lack- 
ing in certain aspects, for example in that there 
were no control field observations. 

The observations were made on the nights 
of 2005, April 18, 22 and 25 with the SIRIUS 
(Simultaneous-three-colour Infrared Imager for 
Unbiased Survey) camera of the Infra Red Survey 
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Facility's 1.4m telescope located at the SAAO's 
Sutherland observatory. SIRIUS has three science- 
grade HAWAII 1024 x 1024 arrays with two di- 
croic mirrors which enable simultaneous observa- 
tions in J (A = 1.25/xm), H (A = 1.65£im) and 
K s (A = 2.25/im). Precise colour information can 
thus be obtained with such simultaneous observa- 
tions. The atmospheric conditions for the three 
nights gave a FWHM for the seeing disk of 1.1 to 
1.4 arcsec for images in the K s band. The field of 
view was 7.8 x 7.8 arcmin 2 with a pixel scale of 
0.45 arcsec per pixel. 

The total exposure of the observations consisted 
of three nights of 12 dithered sets per night, each 
set consisting of 10 ditherings of 30s each. This 
resulted in a total exposure time of 10800 seconds 
for the three nights. The initial data reduction was 
done at the University of Cape Town using the 
SIRIUS pipe-line software which include all the 
basic reduction steps such as dark subtraction, flat 
fielding, and sky subtraction to produce a single 
image for one night's observations. The images of 
the three nights were median combined to produce 
a final image. 

The IRAF task DAOFIND was used to select 
stars from the field. In total 1283 objects were de- 
tected in all three bands. The instrumental magni- 
tudes were calibrated to the 2MASS system. This 
was done by selecting 40 stars in the observed 
field detected in J, H, and K for both the IRSF 
and 2MASS and for which the world coordinates 
of both sets agreed such that the cross identifi- 
cation between the IRSF and 2MASS is unique. 
The correlation between IRSF instrumental and 
2MASS apparent magnitudes as well as between 
IRSF colours and 2MASS colours of the 40 stars 
were investigated. This led to the following equa- 
tions transforming IRSF instrumental magnitudes 
to 2MASS apparent magnitudes: 

J 2 mas S = Jirsf+0.1175(J-H)irsf-- 4.2318 
H 2m a SS = Hirsf - 0.0761 (H - Ks)msF - 4.0974 
Ks 2m ass = Ksirsf + 0.0240(H - Ks)msF - 5.0834 

Following this calibration our observations reached 
J = 19.79, H = 18.43, and K s = 17.64. 

It is also necessary to comment on adopting a 
distance of 2.5 kpc to RCW 34 before present- 
ing the results. I t has already been noted that 
Hevdari-Malaveri ( 19881) asigned a distance of 2.9 
kpc to RCW 34 a lso using the spe ctral type of 
the exciting star. iBik et al. (|2010h pointed out 



that the di fference between their dis tance estimate 
and that ofE evdari-Malaveril (|l9881 ) lies in the dif- 
ferent absolute magnitudes used. Two kinematic 
distance estimates for RCW 34 are also available. 
The water maser associate d with RCW 34 has a 
velocity V kr = 10.8 kms" 1 dBraz k Scaliselll982l) 
Using the rotation curve of IWouterloot fc Brand! 
(|1989l ) this gives a distance of 2.1 kpc. The sec- 
ond kinematic distance comes from the velocity 
of the radio reco mbination line ass o ciated with 
the H II region. ICaswell fc Havned (<1987l) gives 
this as 12 kins -1 which translates to a distance 
of 2.3 kpc, again using t he rotation curve of 
IWouterloot k Brand! (|l989l ). Although a detailed 
analysis of the uncertainties of each of the four 
estimates can be made to decide which to use, it 
seems more appropriate for the present to use the 
mean of the four estimates which is 2.45 kpc. It 
therefore seems r easonable t o sim ply use the dis- 
tance estimate of Bik et al. I (l2010l) . ie. 2.5 kpc. 



3. Results and analysis 

3.1. The two-colour diagram 

The NIR J-H, H-K two colour diagram has 
been used in the past by numerous authors to 
present NIR photometry of young stellar popu- 
lations and to identify different types of objects. 
Here wc follow the same practice. The main re- 
sult of the present study is presented in Fig. [1] 
The main sequence and giant branch were taken 
from Koornneel (Il983h and transformed to the 
2MASS system Four reddening lines usin g the ex- 
tinction law of Rieke k Lebofskv ( 1985! ) are also 
shown. Ri and R2 are the reddening lines that 
bracket the main sequence and giant branch while 
R3 is the reddening line for a star of spectral 
type F0, which is more or less the upper mass 
limit for T Tauri stars. Weak line T Tauri stars 
(WTTS) have intrinsic NIR co lours consistent 
with those of normal dw arf stars (jLada k Adams! 
119921 ; iMever et allll997t) and should therefore fall 
between reddening lines Ri and R 3 . The solid 
green line g ives the class ic al T T auri star (CTTS) 
locus from IMever et al.l (|l997t l and the dashed 
green lines the upper and lower boundaries as cal- 
culated from the errors given by these authors. 
The meaning of reddening line R4 is quite obvi- 
ous. 

We broadly identify three groupings or cluster- 
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Fig. 1. — Near- infrared two-colour diagram for 
1283 objects (red crosses) dete cted in J, H and K 
Blue triangles are WTTS fromlStrom et all (Il989l) 
and black dots CT TS from lStrom et al l(|l989t) and 



Cieza et al.l (|2005l ) 



ings of objects on the two-colour diagram, with 
the area that each more or less cover indicated by 
the ellipses in Fig. [TJ Group A is a strong cluster- 
ing of objects extending along and on both sides of 
reddening line R2 just above its intersection with 
the CTTS locus and is the most obvious aspect 
of the distribution of the 1283 objects on the two- 
colour diagram. Above reddening line R2 most of 
the group A objects lie between reddening lines R2 
and R 3 . There also seems to be a band of objects 
starting from the Ay = 10 point on reddening line 
R 2 and extending almost along the CTTS locus to 
the lower left with most of the objects lying inside 
the boundaries of the CTTS region although some 
lie just above the upper boundary. This is group 
B. Group C is a very loose group of objects that 
lie to the left of group A. The reason for identify- 
ing group C is mainly because there seems to be a 
gap between the bluer objects in group B and the 
objects in group C. 

It is also seen that there are highly reddended 
objects with visual extinctions around 20 as well as 
objects with very large infrared excesses. For the 
present we focus mainly on trying to understand 
group A which is the most obvious feature of the 
distribution on the two colour diagram. 

Given the distribution of objects on the two- 



Fig. 2. — Near-infrared two-colour diagram for 
IRSF objects brighter than J = 16.21 magnitude. 

colour diagram in Fig. [T] the question might be 
asked as to the reality of the clustering and the na- 
ture of the objects belonging to group A. Although 
in NIR studies of other star forming regions it is 
not uncommon to find objects in the same region 
on the two-colour diagram as that of group A, the 
large number of objects in the case of RCW 34, 
is uncommon. The fact that the SIRIUS camera 
has three independent arrays means that the de- 
tection of an object in all three the NIR bands 
implies a real detection and not some artifact on 
the array. The existence of a large number of ob- 
jects with NIR colours as reflected by the group A 
objects should therefore be regarded as real. Since 
we have no spectroscopic identification of any of 
the objects it is necessary to proceed without such 
information. 

Focusing attention on the group A objects, we 
note the following. First, they cluster just above 
the CTTS locus. In fact it would seem as if the 
CTTS locus acts as an approximate lower bound- 
ary for these objects on the two-colour diagram. 
Second, a significant fraction of the objects lie to 
the right of reddening line R2. These objects have 
an infrared excess and cannot be dereddended to 
the main sequence. It is also seen that a further 
significant fraction of the remainder of the group 
A objects lie above the CTTS locus and between 
reddening lines R2 and R3. Although these ob- 
jects can be dereddened to the main sequence, 
such a dereddening would imply the presence of 
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a very large number of stars of spectral type ear- 
lier than FO. Although the visual extinction of 
these stars would be between 10 and 15 magni- 
tudes, the dereddening would also result in some 
being ionizing stars which definitely would have 
been detected as H II regions in radio continuum 
surveys or as bright infrared sources. Only three 
OB stars are, h owever, associated with RCW 34 
ik et alJl2O10h . 

We note that in some studies similar to the 
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presen t one , such as eg. that of iDahm fc Simon 
( 2005 ) and Barentsen et al. ( 201lh . use has been 
made of 2MASS for the NIR photometry. Inspec- 
tion of the magnitudes of the 2MASS objects in 
the same field covered by the present IRSF obser- 
vations showed that the limiting J-band magni- 
tude is 16.12 for objects detected simultaneously 
in J, H, and K. Figure [5] shows the two-colour 
diagram for IRSF detections with J magnitude 
brighter than 16.12. What is important to note 
is the near absence of objects in the group A re- 
gion. On the other hand, it is seen that there 
still are quite a number of objects in group B. 
This suggests that the majority of the group A 
objects is a population of fainter objects revealed 
by our significantly deeper imaging of RCW 34 
compared to 2MASS. As already argued, it seems 
unlikely that the group A objects are simply red- 
dened main-sequence stars. To avoid a completely 
skewed IMF they rather seem to be lower mass 
objects associated with the RCW 34 star forming 
region. 

The fact that the group A objects cluster above 
the CTTS locus is suggestive that they might be 
reddended CTTS. Some early NIR i maging of em- 
bedd ed clusters s uch as eg. IC348 (lLada fc Ladal 



1995), NGC1333 (lLada et al.ll996l) . L1630 dLi et al 



19971 ). as well as some more recent imaging e g. in 
the cases of NGC 2316 dTeixeira et al.ll2004l) and 
the Horschcad Nebula (B owler et al T l2009h do not 
show similar large numbers or even has a lack 
of apparent CTTS's. This may well be due to 
these surveys not being deep enough or that in- 
trinsically there is an absence of large numbers of 
CTTS in these star forming regions. Although not 
in such large numbers as in the case of RCW 34, 
other star forming regions such as eg. NGC 7538 
( Balog et al.l 2004 s ) do show a significant number 
of cluster members which, in terms of Fig. [TJ lie 
just below reddening line R2 and thus can not 
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Fig. 3. — Near-infrared two-colour diagrams for 
the Orion Nebula Cluster (top panel) and IC 348 
(lower panel). 



be dereddened to the main-sequence as is also 
the case for RCW 34. Similarly, for AFGL5180 
( Devine et al.l 2008T ) quite a significant fraction of 
the objects on the two-colour diagram lie below 
th e reddening l i ne R? but above the T Tauri locus 
of iMever et all (|l997f ). 

In all of the above mentioned studies objects ly- 
ing to the right of reddening line R2 and above the 
T Tauri locus are regarded as pre-m ain sequence 
stars. In fact, lLada fc Adamsl (J1992J) have shown 
that pre-main sequence stars classified as CTTS 
based on the equivalent widths of their Ha emis- 
sion lie in the group A regi o n. Other studies eg . 



that of IStrom et al. (Il989h. ICieza et al.l (1200511 1 . 
Luhman et al.l ( 19981 ) , and Barentsen et al. (|2011 ) 
confirm this. In Fig. Q] we also show the po- 
sitions of the spectroscopically identified CTTS 
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from IStrom et all l|l989l ) and ICieza et all (120051) 
as well as the WTTS from lStrom et al.l (jl989h . It 
is seen that the spectroscopically identified CTTS 
are associated with the group A and B regions. 
Luhman et all (|l998f ) also found that for IC 348 
sources showing signs of disk activity are associ- 
ated with our group A objects. 

As examples and for comparison with RCW 34 
we show in Fig. [3] the two-colour diagrams for 
th e Orion Nebula C l uster (ONC), using the data 
of lHillenbrand et al.l (Il998l). and IC 348, using the 
data of lLuhman et al.l (|2003l ) . It is seen that while 
the ONC has quite a number of objects in region 
A (which we will later use), on the other hand IC 
348 has very few. None of the two, however, has 
the same clustering as seen in RCW 34. 

Considering the above, a preliminary conclu- 
sion is that the group A objects are most likely 
lower mass pre-main sequence stars and that some 
of the group B objects, although not having an in- 
frared excess, might also be T Tauri stars. 
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Fig. 4. — Colour-magnitude diagram for RCW 34. 
The red crosses are the objects detected with the 
IRSF. The solid black line is the unreddened main 
sequence for a distance of 2.5 kpc. The dashed 
line is the corresponding main-sequence reddened 
with Ay = 5.1 magnitudes. Blue triangles and 
solid squares are respectively the late-type main 
sequence s tars a nd pre-main sequence stars from 
Bik et all (|2010h . 



3.2. The colour-magnitude diagram. 

In Fig. |4] we show the K vs J-K colour- 
magnitude diagram (CMD) for the 1283 objects 
detected in J. H. and K with the IRSF (red 
crosses). We use the JK CMD in order to compare 
the positions of the spectroscopically identified 
late- ty pe and pre-main sequence stars of lBik et al 



(|2010h (taken from their Table 6) with our NIR de- 
tections. The main-seque nce (solid line) was t aken 
from the Padova models (jMarigo et al.l 120081 ) and 
has been adjusted for distance only using a dis- 
tance of 2.5 kpc to RCW 34. The red dots on 
the main sequence are the positions of the more 
massive main sequence stars with their masses (in 
solar masses) given to the left of each point. We 
also show in Fig. [4] the po sitions of those stars 
for which iBik et aL ( 2010| ) were able to deter- 
mine spectral types from the SINFONI observa- 
tions. The blue triangles are the main sequence 
stars and the solid black squares the late-type pre- 
main sequence stars. Although there is some over- 
lap in the J-K colours of the two groups, it is 
seen that, except for one case, the pre-main se- 
quence stars are on average redder than the main- 
sequence s tars. The reddest object in the list of 



Bik et all (|2010l ) is star 24, with J-K ~ 4. It is 



classif ied as an early K-type dwarf by IBik et al 



(|2010h but of uncertain luminosity class. Its posi- 
tion is also shown in Fig. 21 

It is seen that the objects detected with the 
IRSF populate a very specific area on the CMD. 
Most of the objects have K < 14 and 1 < J - K < 
3 with a clustering around K < 16 and J-K ~ 2. 
The effect of the limited sensitivity of our imag- 
ing is also clearly visible as the sharp cutoff on 
the lower right hand side of the region covered 
by the I RSF objects. Insp ection of the equivalent 
CMD of IBik et all (|2010h (their Figure 7) shows 
that these authors also detected some objects that 
lie between K magnitudes 14 and 16 and for J-K 
just less than 2 there is a clustering of objects to 
the right of the main sequence. Our imaging also 
shows a large number of objects at that position 
on the colour-magnitude diagram. 

Using the v isual extinctions calculated by 
Bik et al. | (|2010h (their Table 4) we also calculated 
the dereddened positions of the three early-type 
stars as well as that of star 24. The dereddening 
vectors are shown by the arrows starting at each 
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of the stars. For the three early-type stars the 
dereddended positions do not fall exactly on the 
main sequence but slightly to the right. Adding a 
further approximately 0.9 visual magnitudes will 
bring the three dereddened early-type stars on the 
main sequence. However, we noted that there is 
a difference in extinction for t he three OB star s 
as given in Tables 4 and 6 of iBik et all (|2010f ). 
For example for star 1, the exciting star, Tabic 
4 gives an extinction of 4.2 while it is given as 
4.8 in Table 6. The reason for this not clear. We 
also note t hat th e extinction toward region II of 



Bik et al 



(see their Figure 4), which is 
close to the location of the three OB stars, is A„ 
= 5.1 magnitudes. The required visual extinction 
of 5.1 magnitudes to bring the unreddened main 
sequence to the positions of the three OB stars 
is therefore i n general agreem ent with the mea- 



surements of Bik ct al 



( 2010l ). It is furthermore 
interesting to note that star 24, with an estimated 
extinction of 15.9 ± 2.8 cannot be dereddened at 
all to the main sequence and if it could, it would 
not be a K-dwarf. Star 24 is therefore most likely 
still in the pre-main sequence phase. 

Taking the three OB stars' positions to indicate 
the position of the reddened main sequnce, we can 
redden the unreddened main sequence accordingly. 
This is shown as the dashed line in Fig. |4] It is 
seen that the majority of our IRSF sources still 
lie significantly to the right of the main sequence 
as will be expected for pre-main sequence stars. 
This result is in support of our hypothesis that 
the sample of 1283 objects detected in J, H, and K 
contains a significant number of pre-main sequence 
stars. Those objects lying significantly to the left 
of the reddened main sequence are most probably 
foreground stars. 

3.3. The luminosity function 

The luminosity function is an important char- 
acteristic of the population of pre-main sequence 
stars since it contains information about the age, 
initial mass function, and star formation history of 
that particular population. In constructing the lu- 
minosity function and estimating the ages of mem- 
bers of the population the ideal would be to know 
the effective temperature and luminosity of each 
individual star from spectroscopic measurements 
and to then construct a luminosity-temperature 
diagram on which model isochrones can overlaid 



to estimate the pre-main sequence ages. In the 
present case we only have NIR photometry data 
which is not suitable to determine stellar effective 
temperatures due to contamination from excess 
emission from circumstellar material. Uncertain- 
ties in the photometry can also lead to erroneous 
conclus ions about the a ges of individual objects 
(see eg iPreibischl (|2012h for an extensive discus- 
sion) . 

However, as argued bv lBontemps et ah ( 2001 ). 
the J-band flux is least affected by contamina- 
tion from circumstellar emission and is close to the 
peak of the photospheric spectral energy distribu- 
tion for these cool stars and may be used to es- 
timate stellar luminosities. We therefore followed 
these authors and estimated the stellar luminos- 
ity directly from t he absolute J-band m agnitude 
using the relation (jBontemps et al]|2001l ) 



log 10 (£*) = 1.49 -0.466 Mj (1) 

The estimated uncertainty on log 10 (i*) using Eq. 
□ is 0.19 dex. 
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Fig. 5. — Bolometric luminosity distribution for 
a subset 745 infra-red excess sources (black dots). 
The blue, dashed black, and red lines are the sin- 
gle stellar population luminosity functions respec- 
tively at 10 6 and 10 7 years and for the ZAMS as 
calculated from the mass-luminosity relations at 
these times as given by the models of ISiess et al. 
(<2000h . The luminosity function for Orion is given 
by the black line with red dots. 

Application of Eq. [1] requires calculation of the 
absolute J-band magnitude and therefore dered- 
dening of the observed J-band magnitudes. The 
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assumption was therefore made that all the group 
A objects are reddened CTTS. However, we se- 
lected only the subset of 745 group A objects that 
lie to the right of reddening line R3, above the T- 
Tauri locus and with 0.53 < H - K < 1. These 
objects all have an infrared excess and cannot be 
dereddened to the main sequence. Since we do not 
know the true unreddened intrinsic colours, each 
object was dereddened to a random position be- 
tween the upper and lower boundaries of the T 
Tauri locus but still lying on its individual dered- 
dening line. One instance of dereddening all 745 
objects then results in a luminosity function for 
the given random positions around the T Tauri 
locus. By repeating the procedure a large number 
of times it is possible to construct an average lu- 
minosity function that is representative of a large 
number of realizations of intrinsic colours for the 
dereddened objects. All objects were assumed to 
lie at a distance of 2.5 kpc to calculate the absolute 
magnitudes and therefore the stellar luminosities. 

In Fig. [5] we show with the black solid line 
the luminosity function averaged over 1000 dered- 
dening realizations and normalized such that 
J (dN/dL)dL = 1. Normalization is necessary 
to be able to compare the luminosity functions 
of different star forming regions as we will do be- 
low. The distribution is seen to have a peak at 
about 0.04L Q below which there is a sharp cut-off 
which most likely is a sensitivity effect. Except 
for the last three highest luminosities, the lumi- 
nosity function follows a power law with an index 
of —2.34 ± 0.07 for luminosities above about 0.3 
L . 

The observed luminosity distribution is a re- 
sult of the combination of the evolution of the 
pre-main sequence mass-luminosity relation with 
time, the underlying initial mass function and 
the star formation history of the region. We 
attempted to interpret the observed luminosity 
distribution with single stellar population (SSP) 
luminosity functions constructed by using the 
mass-luminosity relation for pre-main sequence 
isochron es at 10 6 and 10 7 years and the ZAMS 
from the Siess et al. (2000) models combined with 
a Kr oupa stellar initial mass function ( Kroupal 
20011 ) for masses between 0.3 M Q and 2 M . 
For an SSP the luminosity function, dN/dL, is 
related to the initial mass function, dN/dm, as 
dN/dL = (dN/dm) x (dm/dL). To calculate 



dm/dL we approximate d the mass- l umin osity re- 
lation, as given by the ISiess et al model 
for a specific isochrone, on a log-log scale with a 
polinomial which is easily differentiable. The use 
of the 10 6 and 10 7 year luminosity functions as 
comparison with the observed luminosity function 
should not be interpreted as attemps for abso- 
lute age estimates but as guides to interpret the 
observed luminosity function. The result of this 
calculation is shown in Fig. [5]where the luminosity 
functions at 10 6 (blue solid line) and 10 7 (dashed 
black line) years have been adjusted vertically to 
more or less coincide with parts of the observed 
luminosity function of RCW 34. The luminosity 
function for the ZAMS (solid red line) has been 
placed at an arbitrary position. 

Inspection of Fig. [5] shows that although the 
observed luminosity function can be fitted by a 
single power law for L > 0.3 Lq, neither of the 10 6 
or 10 7 year SSP luminosity functions can explain 
the entire observed luminosity function. It is seen 
that the theoretical 10 6 year luminosity function 
explains the observed luminosity function rather 
well for luminosities greater than about 1.4 Lq. 
Below ~ 1.4 Lq the theoretical luminosity dips be- 
low the observed luminosity function, suggesting 
the presence of an older component. The theoret- 
ical luminosity function for an SSP with an age 
of 10 7 years is seen to follow the observed lumi- 
nosity function quite well between ~ 0.2 Lq and 
just less than 1.4 Lq. For L < 0.2 Lq the observed 
luminosity function starts to bend away from and 
fall below the theoretical 10 7 year luminosity func- 
tion. This behaviour is most certainly due to the 
limited sensitivity of our imaging. Finally we note 
that the observed luminosity function nowhere has 
the behaviour of the luminosity function for the 
ZAMS suggesting that the group A objects used 
here are not main sequence objects. 

The fact that a significant part of the RCW 34 
luminosity function can be explained by a com- 
bination of the 10 6 and 10 7 year SSP luminos- 
ity functions is suggestive that star formation in 
RCW 34 has been an ongoing process for about 
10 7 years. As already noted, based on the pres- 
ence of the Q8.5V star that powers the H II 
region, Bik et al. (|2010l ) put an upper limit of 
8 x 10 6 ~ 10 7 years on the age of RCW 34. The 
presence of a class II methanol maser indicates a 
much more recent episode of star formation ac- 
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tivity associated with RCW 34 since these masers 
are known to be exclusively associated w ith a very 
early phase of massive star formation (jEllingsenl 
2006). Furthermore, the 6.7 GHz methanol maser 
lifetime is estim ated to be between 2.5 x 10 4 and 
4.0 x 10 4 years (jvan der Waltl 120051 ) . From their 



3.4. Spatial distribution of objects. 



set of 26 spec troscopica ll y iden tified pre-main se- 
quence stars Bik et al. ( 2010l ) estimate the age 
of the youngest pre-main sequence objects to be 
(2 ± 1) x 10 6 years. Our analysis of the NIR pho- 
tometric data of a completely different sample of 
objects in RCW 34 is there f ore in general agree- 
ment with that of iBik et all (|2010l) . 

In Fig. [5]we also show the luminosity function 
for 204 pre-main sequence star s from the Orion 
Nebul ar Cluster using the data of Hillenbrand et all 
(|l998l) . Th e JHK colours of a l l the objects in the 
catalog of iHillenbrand et al. (1998) were trans- 
formed from the CIT to the 2MASS system and 
afterwards subjected to the same selection criteria 
that was used to select the subset of 745 group A 
sources from which the RCW 34 luminosity func- 
tion constructed. Dereddening was done in the 
same way as for the RCW 34 sources and an aver- 
age luminosity function was calculated over 1000 
dereddening realizations. 

For the Orion sources the luminosity function 
peaks at about 0.1 L Q compared to 0.04 L Q for 
RCW 34. The Orion luminosity function crosses 
the RCW 34 luminosity function between 0.2 and 
0.3 L Q . For L > 2 Lq the two luminosity func- 
tions runs almost parallel to each other except for 
the last two points for RCW 34. In fact, the main 
trend of the ONC luminosity function for L > 2 Lq 
is a power law with index -2.17 ± 0.05 which is 
very similar to that of RCW 34. It is also seen that 
for L > 2 L Q the 10 6 year SSP luminosity function 
(dashed blue line) describes the Orion luminosity 
function quite well. This is in agreement with in- 
dependent age estim ates of about 2 x 10 6 years 
for the ONC (see eg. Reggiani et al. 1 2 1 lh . Com- 
parison of the 10 6 and 10 7 year SSP luminosity 
functions suggest that the slope of the luminosity 
function for L > 2 L© is definitely dependent on 
the age of the system. Thus, just by comparison 
of the luminosity function of RCW 34 with that 
of the ONC already suggests the existence of com- 
ponent with an age of about 1 - 2xl0 6 years in 
RCW 34. 
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Fig. 6. — Top panel: The spatial distribution of all 
1283 objects. Bottom panel: Spatial distribution 
for 91 objects with L > 0.5 L Q (blue dots) and 209 
objects with L < 0.075 L (red dots). The black 
dot indicates the position of the 6.7 GHz methanol 
maser and the star the position of the 08.5V star. 
The circle in the top and bottom panels indicates 
the approximate size and location of the opacity 
hole (bubble) as revealed in the dust emission o f 
the region (See eg. Figs. 1 & 3 of lBik et alJl20ld ). 
Note that the coordinates are in decimal degree. 

In the top panel of Fig. [S] we show the spa- 
tial distribution of the 1283 objects detected in J, 
H and K. Some degree of clustering can be seen 
in a short band running north-east to south-west 
just south of the central 08.5V star. Otherwise, 
the detected objects seems to be more or less uni- 
formly distributed over the field. 

Given that our analysis of the luminosity func- 
tion suggests that star formation in RCW 34 took 
place over an extensive period of time, the ques- 
tion is whether there is any difference in the spa- 
tial distribution of the younger and older pre-main 
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not corrected the larger annulli for the smaller sur- 
face area enclosed inside the frame borders. 



Reggiani et al.1 (|201ll ) recently also investigated 
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Fig. 7. — Radial surface density distribution of 
group A objects with luminosity greater than 
0.5Lq (red line) and objects with luminosity less 
than O.O75L (black line). The origin was taken 
as the position of the 08.5V star. 



sequence stars. Here we follow the discussion in 
the previous section and use luminosity as an ap- 
proximate indicator of age with younger objects 
having greater luminosity than older ones. To en- 
sure that we have a statistically sufficient number 
of objects we select the brighter objects as those 
with L > 0.5 L Q . This resulted in 91 objects. For 
the older group we selected the 209 objects with 
L < 0.075 Lq. The spatial distribution of the two 
groups is shown in the bottom panel of Fig. [6] 

Inspection of the bottom panel of Fig. [5] sug- 
gests a clustering of the higher luminosity objects 
(younger ones; blue dots) closer to the 08.5V star 
than the fainter objects (red dots) and in particu- 
lar in the region of the opacity hole. To quantify 
the distributions we calculated the relative surface 
density of younger and older objects as a function 
of distance from the 08.5V star using a number of 
annulli. The two distributions are shown in Fig. 
[7J The brighter objects clearly shows a signifi- 
cantly higher probability of being found closer to 
the 08.5V star than the fainter objects which have 
a flatter surface density distribution. In both cases 
the decrease in the surface density beyond ~240 
arcsec is artificial and due to the fact that we have 



the variation of the surface density distribution 
of pre-main sequence stars with age in the ONC. 
Interestingly enough these authors found that the 
older objects have a less concentrated distribution 
compared to the younger objects. Qualitatively 
the same behaviour is seen in RCW 34. However, 
it should be kept in mind that the older lower mass 
stars are most probably of the same age as the 
three OB stars. The concentration of the younger 
group closer to the three OB stars seems to suggest 
that the more recent episode of star formation may 
have been triggered by the presence of these three 
stars. 

4. Summary and Conclusions 

We presented NIR imaging data on RCW 34 
for a 7.8 x 7.8 arcmin 2 region centered on the 6.7 
GHz methanol maser associated with RCW 34. A 
total of 1283 objects were detected in J, H, and 
K-bands. The distribution of these objects on the 
two-colour diagram shows a concentration of more 
than 750 objects for which the colours are the same 
as that of confirmed classical T Tauri stars found 
in other star forming regions. Given that the po- 
sition of the main sequence on the K vs J-K CMD 
is determined by the positions of the three OB 
stars, the distribution of the IRSF sources on the 
CMD is also suggestive that a significant number 
of these are lower mass pre-main sequence stars. 

We also constructed the bolometric luminos- 
ity function for the 745 objects and showed that 
parts of the luminosity function can be explained 
by SSP luminosity functions with ages of 10 6 and 
10 7 years. The presence of a young component 
in the stellar populatio n of RCW 34 is in agree- 
ment with the results of Bik et al. ( 2010l ) based on 
26 spcctroscopically identified pre-main sequence 
stars as well as the presence of a class II methanol 
maser which indicate a very recent episode of mas- 
sive star formation. Our estimate of 10 7 years for 
the age of an older component is qualitatively in 
agreement with the main sequence life times of the 
three OB stars associated with RCW 34. 

Whereas previous studies of RCW 34 focused 
more on the central region around the ionizing 
star, our NIR imaging revealed a more spread out 
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population of low mass pre-main sequence stars. 
The younger stars appear to be more concentrated 
in the central region around the 08.5V star while 
the older pre-main sequence stars seems to be 
more spread out. 

Consid e ring o ur NIR results as well as that of 
Bik et all (|2010h . RCW 34 seems to be a much 
more interesting star forming region than was per- 
haps previously thought. First, there appears to 
be a very large number of low mass stars formed 
over a period of about 10 7 years with the older 
being more spread out than the younger compo- 
nent. Second, it seems rather certain that mas- 
sive star formation in RCW 34 took place in two 
separate events with only three OB stars form- 
ing in the first event and most probably only a 
single massive star, as evidenced by the 6.7 GHz 
maser, forming in the most recent event. The fact 
that the brighter (younger) lower mass pre-main 
sequence stars seems to cluster around the posi- 
tion of the three OB stars strongly suggests that 
that the physical effects these three stars had on 
the remainder of the molecular cloud from which 
they formed could have played a role in the more 
recent episode of star formation. Obviously our 
photometric study need to be followed up by a 
spectroscopic study which we plan to do. Apart 
from verifying our photometric identifications, will 
such a study be very usefull to also try to unravel 
the star formation history, especially of the low 
mass stars, in RCW 34. Given that there also 
seems to be a significant number of low mass pre- 
main sequence stars spread out almost uniformly 
over the 7.8 x 7.8 arcmin 2 region might require a 
different scenario for the star for mation history i n 
RCW 34 than that suggested bv lBik et all (|2010h . 
Although RCW 34 cannot be regarded as peculiar 
anymore, it is seems to be interesting and suffi- 
ciently different from other star forming regions 
to require further investigation. 

We would like to thank an anonomous referee 
for constructive comments to improve the paper. 
This work was supported by the National Research 
Foundation under Grant number 2053475. 
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